Density Matrix Renormalization Group Study of Incompressible Fractional Quantum Hall States 
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We develop the Density Matrix Renormalization Group (DMRG) technique for numerically studying incom- 
pressible fractional quantum Hall (FQH) states on the sphere. We calculate accurate estimates for ground state 
energies and excitation gaps at FQH filling fractions v = 1/3 and v = 5/2 for systems that are almost twice as 
large as the largest ever studied by exact diagonalization. We establish, by carefully comparing with existing nu- 
merical results on smaller systems, that DMRG is a highly effective numerical tool for studying incompressible 
FQH states. 
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Introduction - Numerical techniques, particularly the ex- 
act diagonalization of a small system of N e interacting 2D 
electrons in a strong perpendicular magnetic field, have been 
a central and essential tool in the theoretical study of the frac- 
tional quantum Hall effect (FQHE) [1] dating back, rather cu- 
riously, to the days |2j,l3[] before Laughlin wrote down [4] his 
famous wavefunction. This wavefunction has served as the 
paradigm for understanding the odd-denominator incompress- 
ible FQH states during the last 25 years due, in part, to sup- 
porting evidence given by exact diagonalization studies: the 
Laughlin wavefunction has almost unity overlap with the nu- 
merically calculated exact many-body wavefunction for small 
(N e = 4 — 8) spherical systems interacting via the Coulomb 
interaction, as shown in the seminal work by Haldane H. Of 
equal importance, Haldane developed an extremely enlight- 
ening pseudopotential expansion in the conserved relative or- 
bital angular momentum channel, establishing in the process 
the amazing fact that the Laughlin wavefunction is the exact 
eigenstate for the short-ranged interaction potential V\ ^ 0, 
Vj = for J = 3,5,..., where the Vj is the interaction be- 
tween two electrons with relative angular momentum J. This 
immediately led to a deep physical and mathematical under- 
standing of why and when a Laughlin-type incompressible 
state would be a good description of the ground state observed 
in experiments. Ever since Haldane's work, exact diagonal- 
ization of finite systems, mostly on the sphere, torus, or disc, 
has been a crucial theoretical tool for studying FQH systems - 
indeed, the only one with great predictive power for the depen- 
dence of the ground state on material parameters. It is there- 
fore of substantial interest and importance that new numerical 
methods are developed for larger system sizes which could be 
useful in understanding more complex and exotic FQH states. 
In particular, there is great current interest in higher Landau 
level (LL) FQH states in the first |0] or even the second 0] 
LL and the fascinating even denominator states (e.g. the 5/2 
FQH state [8] discovered 20 years ago) as well as studies of 
spin polarization properties |9[] of FQH states. These all re- 
quire finite size studies involving system sizes N e much larger 
than what has been achieved through the exact diagonalization 



method, which is typically restricted to N e = 10 — 18 depend- 
ing on the LL filling factor v. Since the Hilbert space dimen- 
sion grows exponentially with N e , the exact diagonalization 
technique has essentially reached its (memory and storage) 
limit, and no significant advance in increasing N e for exact 
diagonalization is anticipated in the near future. 

In this Letter we use a powerful technique, the Density Ma- 
trix Renormalization Group (DMRG) |T3 HI, for studying 
ground and excited state properties of incompressible FQH 
states, allowing a substantial jump (~ by almost a factor of 
2!) in N e (up to N e = 20 for v = 1/3 and 30 for v = 5/2). 
We introduce, establish, and validate this technique by con- 
centrating on the archetypal FQH state at v = 1/3 and the 
mysterious and not yet well-understood FQH state at v = 5/2. 
Our work not only reproduces the numerical results obtained 
earlier in exact diagonalization studies for the same values of 
N e , but also extends our understanding of the ground state 
and low-lying excited state properties of these FQH states to 
substantially larger system sizes. Our treatment of large sys- 
tems on the sphere depends on the use of intricate, non-trivial 
methods which are described below, unlike previous pioneer- 
ing applications of the DMRG to the quantum Hall problem, 
which concentrated mostly on compressible non-quantized 
states (e.g. bubbles and stripes) in higher Landau levels on 
the torus 11211 . We suggest that these techniques will allow, in 
the future, DMRG studies of problems simply impossible to 
treat with exact diagonalization. We believe that the develop- 
ment of the DMRG method (which is almost exclusively used 
for ID interacting systems) to the 2D FQHE problem should 
have high impact on both the physics of the FQHE and the 
theory of the DMRG. 

DMRG - The Density Matrix Renormalization Group 
(DMRG) EEl method is a sophisticated algorithm to cal- 
culate the ground-state and low-lying excitations of quantum 
many-body problems. It can treat very large systems with 
hundreds of spins or electrons, and it is extremely accurate 
for quasi one-dimensional systems. The basic idea consists of 
truncating the Hilbert space in such a way as to maximize the 
amount of information we preserve about the actual ground 
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FIG. 1 : a) Ground state and variational energies at v — 1/3 obtained 
with DMRG and exact diagonalization (ED). The extrapolated values 
are obtained using a quadratic function of 1 /N e . Variational results 
are barely visible because they lie on top of the DMRG energies at 
this scale, b) Finite size scaling using a renormalized magnetic length 
In. following Refs. (2fll2lll . 



state of the system. In order to do this, we need to systemati- 
cally apply a change of basis called the density matrix projec- 
tion. The resulting wavefunction is a very accurate representa- 
tion of the ground-state in a reduced Hilbert space with a given 
number of states, and assumes the form of a matrix product 
state. In this paper we want to apply the DMRG to study a 
two-dimensional system of electrons in a strong a magnetic 
field. The generic form of the Hamiltonian is 



H =~ ( m l' m 2l^l TO 3,™4)cL 1 cL 2 Cm 4 Cr, 



(1) 



mi ? rn 2 

m.3 ,7714 



where the indices m represent the quantum numbers of free 
electron orbitals in a given Landau level (LL), and V is 
the electron-electron interaction (which will simply be the 
Coulomb repulsion in this paper although one could modify 
the interaction in order to account for finite layer thickness, 
for instance). Energies are quoted in units of e 2 /Airelo, where 
la = (Hc/eB) 1 / 2 is the magnetic length. The orbitals on a 
spherical geometry are labeled by their angular momentum, 
while on the torus, one needs to use the magnetic translations, 
identified by the two components k x , and k y . We have omitted 
the spin degrees of freedom due to the large Zeeman splitting. 
Thus, this Hamiltonian represents a two-dimensional system 
of spinless electrons with long-range interactions. In order to 
apply DMRG, one can use a generalization of the algorithm 
to momentum space, which has already been applied to the 
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FIG 2: Excitation energy at v = 1/3 obtained with DMRG and 
exact diagonalization, and the extrapolation using a linear function 
of 1/N e . We also show the results for the rescaled energies using a 
size-independent magnetic length. 



Hubbard model|[l3|. We note that the extension of the DMRG 
to 2D systems is an important challenge which is a subject of 
great current interest [11]. Our method exploits the Landau 
level structure of the FQH problem by applying the momen- 
tum space formulation of the DMRG. 

The idea is to arrange the orbitals along a one-dimensional 
chain, and apply the finite-system version of the DMRG. The 
simulation starts with a small system with a few orbitals that 
can be easily diagonalized. We split the system in two halves: 
the left block (Bl) and the right block {Br). We increase the 
system size by systematically adding orbitals to the system. 
We do that by including two orbitals in the center, between the 
blocks, in a configuration B^ • We then calculate the 

ground state of the Hamiltonian, and obtain the reduced den- 
sity matrix of one block plus the new added site, say [£>£•], by 
tracing over the basis states of [•-£?«]. We pick the new basis 
for the new block Bl as the eigenstates of the reduced density 
matrix with the largest eigenvalues. We can adjust how many 
basis states we keep, thereby controlling the accuracy of the 
computation. Once the desired system size has been reached, 
we grow one block and shrink the other one, sweeping from 
left to right, and right to left, until we converge. The mo- 
mentum space version of the DMRG allows us to use the total 
momentum (angular momentum on the sphere), besides the 
particle number N e , to work in a subspace of reduced dimen- 
sion. The matrix elements of the interaction V on the sphere 
take the form |H [3] 

(■mi,m 2 \V\m3,m4) = 

2S +.J 

Y (S,m 1 ,S,m 2 \JM)(JM\S,m 3 ,S,m A )V J /R, 



J=0 M=-J 



(2) 



where R is the radius of the sphere, and the coefficients Vj 
are Haldane's pseudopotentials [5], derived from the pure 
Coulomb potential. They represent the interaction energy of 
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FIG. 3: a) Ground state energies at v = 5/2 for DMRG and exact 
diagonalization. b) Finite size scaling using a renormalized magnetic 
length £ , following Refs.H^IH. 



two electrons with relative angular momentum J, and depend 
on the Landau level. To obtain this expression, we have used 
addition of angular momenta, and 25 = N<$> is the total flux 
through the system (not to be confused with electron spins, 
which are assumed to be fully polarized), given by the num- 
ber of electrons N e and the filling fraction v. As shown by 
Haldane fl, the Laughlin state at v = 1/3 is the exact ground 
state for a potential obtained by setting all Vj = except for 
V\. We use this observation to compute the Laughlin state in 
the same basis as the ground state of the system for Coulomb 
interactions. In order to represent both states in the same ba- 
sis simultaneously with equivalent accuracy we calculate the 
ground states of the full Coulomb interaction, if>, and Hal- 
dane's Hamiltonian, ^n, and we use both states to obtain the 
reduced density matrix at each DMRG step [10]. This is es- 
sential for us to estimate variational energies as well as over- 
laps between the Laughlin state and the exact ground state of 
the system. 

Optimization strategies - Since the Laughlin and Moore- 
Read states have non-zero amplitude only on "squeezable" 
configurations I1 1611 . we have replaced the traditional DMRG 
warmup scheme by starting our simulation from a small set 
of "squeezable" states, thereby improving convergence con- 
siderably. This also removes undesired effects, such as stick- 
ing to an excited state. Another important ingredient of our 
algorithm is the addition of density matrix corrections 01711 
that takes into account the incompleteness of the environment 
block and the missing entanglement due to the non-locality of 
the interactions. These corrections improve the energies by 
one order of magnitude compared to the conventional algo- 
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TABLE I: Ground state energies E1/3 and variational energies Eo 
for v — 1 /3 obtained with DMRG. The last column shows the over- 
lap between the exact wave function ip and the variational (Laughlin) 
state ipo obtained using Haldane's Hamiltonian. 



rithm. These optimizations have allowed us to obtain ground 
state energies and gaps for systems with up to N e = 20 elec- 
trons at v = 1/3 and N e = 26 at v — 5/2, using up to 5000 
DMRG states, with a typical truncation error of the order of 
10~ 6 for the largest systems studied here. It is interesting to 
note that the size of the DMRG Hilbert space is of the order 
of 1, 000, 000 states, while the full Hilbert space has a total 
dimension of the order of many millions for N e > 12. 

Results - In Table I, we show the variational energy of 
the Laughlin state compared to the actual ground state en- 
ergy of the full Coulomb Hamiltonian at v = 1/3, at flux 
N<$, = 3(N — 1). These were obtained using exact diago- 
nalization for N < 12, lfl8ll and DMRG for larger systems. 
We also show the overlap between both wavefunctions, in- 
dicating unambiguously that the Laughlin state does, indeed, 
describe the physics of the v = 1/3 FQHE as we approach the 
thermodynamic limit. The energies are also plotted in FigJT] 
Extrapolations to the thermodynamic limit can be carried out 
using a second-order polynomial in 1/N e , see Fig. 01 a X or a 
first-order one after renormalizing the magnetic length £0 to 
take into account the curvature of the sphere, which puts all 
the data points on a straight line lfl9l, I20T 21 ] . Notice that the 
scale shown in the y-axis in FigQJb) is much smaller than the 
original, making the DMRG errors visible for large systems. 
However these errors are small, of the order of 10~ 5 , and one 
order of magnitude smaller for the variational energy of the 
Laughlin state. 12211 It can clearly be seen that finite-size ef- 
fects are very small for N > 10. The fact that accuracy and 
convergence are better for the Laughlin wavefunction can be 
attributed to the short range of the interactions and the fact 
that it spans a reduced Hilbert space [ 16], allowing DMRG to 
efficiently represent it as a matrix product state 12311 . The ex- 
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/3 



-0.41010 



for the Laughlin state and exact ground state energies, re- 
spectively, agree well with previous Monte Carlo 11911 and 
exact diagonalization|20] estimates of Eq = —0.40975, and 
E1/3 = -0.4101. 

In Fig|2] we show DMRG results for the energy gap at 
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FIG. 4: Size dependence of the excitation energies and gap at v 
5/2, obtained with DMRG. 



v = 1/3, calculated from the energy of charged excitations 
with flux N$ = 3(N — 1) ± 1. Our extrapolated value in the 
thermodynamic limit of A1/3 = 0.1010 agrees well with pre- 
vious estimates of A1/3 = 0.1012 using smaller systems up 
to N e = 12. HH] 

Having successfully tested the DMRG at v = 1/3, we pro- 
ceed by extending the analysis to the second Landau level, at 
v = 5/2. In Fig|3ja) and (b), we show the ground state en- 
ergies at flux N$ = 2N — 3, and the scaled corrected values 
using a size-independent magnetic length, respectively. Ex- 
trapolations yield a value of E 5 / 2 = —0.364, which compares 
well with Morf's result B20I1 of —0.366. Notice that in our 
study we have ignored system sizes smaller that N e = 10 be- 
cause of strong finite-size effects. 

The calculation of the gap at v = 5/2 presents an extra 
complication due to 'aliasing' with states corresponding to 
conventional FQH fractions B20l,l2lll . This constrains the num- 
ber of accessible states to the sequence N e — 10, 14, 18, .... 
Due to numerical limitations, we have been able to add a sin- 
gle extra point to this sequence with N e = 22. The results are 
depicted in Fig0J where we show the energies for the half- 
flux-quantum (e/4) quasiparticle and quasihole, in addition 
to the gap, which is their average. As may be seen from the 
figure there is more scatter in the quasihole and quasiparti- 
cle energies than in the gap. In particular, the energies at 
N e = 18 look anomalous. This suggests strong finites-size 
effects, and an alternative extrapolation based on only the two 
largest system sizes would give a significantly larger extrap- 
olated gap value. Details will be presented elsewhere. The 



results for smaller system sizes coincide with Morf's [20, 21], 
and a linear extrapolation using all the datapoints yields a 
value of the gap of A 5 / 2 = 0.030, also larger than the result of 
A5/2 = 0.025 previously reported. We emphasize that our es- 
timated gap is an upper bound to the real gap since, following 
Morf l20il2in . we treat the quasiholes as point particles. 

Conclusion - We have developed and validated the 
DMRG technique as a theoretical tool for the numerical study 
of FQH incompressible states, thereby increasing the possi- 



ble tractable system size N e by a nearly a factor 2 (e.g. from 
N e = 10 to N e = 20 electrons for v = 1/3). We have 
applied the technique to v = 1/3 and 5/2, obtaining in the 
process the best estimates for both the ground state energy 
and the excitation gap for these important FQH states. We 
have also significantly expanded the regime of validity for the 
infinite system extrapolation, 1/N e — > 0, which is essential 
for any comparison with experimental work. We have found 
that a reliable finite size scaling requires a large number of 
DMRG states (up to 5000 in our work), and that density ma- 
trix corrections and "squeezing" are essential to improve the 
accuracy and convergence. Furthermore, this technique can 
be used to calculate and validate variational states, such as the 
Laughlin state. Our work establishes the DMRG as a numeri- 
cal tool for future theoretical FQHE studies where large values 
of N e are essential for theoretical understanding. In particu- 
lar, the DMRG technique may be able to to shed light on var- 
ious mysterious fractions such as 5/11, 6/13, 4/13, 3/8, and 
should be useful in studying incompressible states in higher 
LLs, where large system sizes are essential in distinguishing 
between competing compressible and incompressible phases 
1 24J] which lie close in energy. 
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